#
# 
#
# Authors: Soeren Gebbert, soeren.gebbert@vti.bund.de
#          Axel Don, axel.don@vti.bund.de
#
# Copyright:
#
# Johann Heinrich von Thünen-Institut
# Institut für Agrarrelevante Klimaforschung
#
# Phone: +49 (0)531 596 2601
#
# Fax:+49 (0)531 596 2699
#
# Mail: ak@vti.bund.de
#
# Bundesallee 50
# 38116 Braunschweig
# Germany
#
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; version 2 of the License.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#

# Set the filename of the output graphical file
pdf("scamp_run.pdf")

# Read the data
run_data  = read.csv("run.csv")

# Set graphical layout
par(mfrow=c(4,1))

# Plot data
plot(micr_mass~timestep, data=run_data, main="Mass of micro-organism over time steps", type="l")
plot(eaten~timestep, data=run_data, main="Consumed carbon rate per time step", type="l")
plot(carb_total_mass~timestep, data=run_data, main="Total carbon particles mass over time steps", type="l")
plot(carb_num~timestep, data=run_data, main="Total number of carbon particles over time steps", type="l")

# Prepare the device
dev.set(dev.next())
